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Consider a stabilizer state on n qudits, each of dimension D with D being a prime or a squarefree 
integer, divided into three mutually disjoint sets or parts. Generalizing a result of Bravyi et al. [J. 
Math. Phys. 47, 062106 (2006)] for qubits [D = 2), we show that up to local unitaries on the three 
parts the state can be written as a tensor product of unentangled single-qudit states, maximally 
entangled EPR pairs, and tripartite GHZ states. We employ this result to obtain a complete 
characterization of the properties of a class of channels associated with stabilizer error-correcting 
codes, along with their complementary channels. 
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I. INTRODUCTION 

The study of entangled quantum states of systems con- 
sisting of two or more parts is a central problem in quan- 
tum information theory. The Schmidt decomposition 
provides a fairly complete characterization of the pure 
states of a bipartite system. However, mixed states on 
bipartite systems and pure states on systems of three or 
more parts present a much more difficult problem — see 
yl for a comprehensive review — and a relatively complete 
understanding of the situation exists only for some very 
special cases. 

The present paper considers the special case of (pure) 
stabilizer states on n qudits, each of dimension £), and 
addresses the problem of characterizing the correspond- 
ing tripartite state when the n qudits are partitioned into 
three disjoint sets A, B, and C, and arbitrary unitary 
transformations are allowed on each of the three parts. 
The case of qubits, D = 2, was studied by Bravyi et al. 
[3], who showed that such a stabilizer state is equivalent, 
up to local unitaries on the three parts, to a tensor prod- 
uct of pure unentangled single qubit states; maximally 
entangled two-qubit states or EPR pairs, with one qubit 
in one part and the other qubit in a different part; and 
GHZ states on three qubits, one lying in each part. In 
this paper we generalize these results to the case D > 2, 
where D is either a prime or a squarefree integer (i.e., not 
divisible by the square of any integer greater than 1). 

The stabilizer formalism 0-l3 was first introduced to 
simplify the construction and analysis of quantum error 
correction codes. Soon thereafter it was generalized from 
qubits to higher dimensional qudits 6, 7!| . Most of the 
codes known when the formalism was introduced, and the 
majority of those discovered since, are stabilizer codes. 
The formalism has also been used for measurement-based 
quantum computation [8| and fault-tolerant topological 
quantum computation [9|. There has been a lot of re- 
search on singie-qubit local unitary (LU) and single- 
qubit local Clifford (LC) e quivalence of qubit stabilizer 
states/graph states pi]| - [l2| but here we consider parti- 
tionings where each part can have several qudits and ar- 



bitrary gates acting on qudits belonging to the same part 
are permitted. 

In [13| we studied a class of channels obtained from qu- 
dit stabilizer (equivalently, additive graph) codes where 
a subset of the carrier qudits is lost. We fully charac- 
terized their information carrying capacities in terms of 
subset information groups, a concept related to the no- 
tion of correctable algebras introduced in [1J|. We also 
provided an efficient algorithm to find the subset infor- 
mation group. In this paper we adopt the name stabilizer 
code channels for such channels. 

The paper is organized as follows: Section |lT] intro- 
duces various concepts that will be used later: Pauli and 
Clifford operators, one- and two-qudit gates, stabilizer 
and graph states. It also contains some mathematical 
results, one of which. Corollary [SJ is of some interest by 
itself: it allows the decomposition of stabilizer states into 
a tensor product of such states when D — did2 • • • is a 
product of mutually coprime factors. In the following 
Section [ml we prove that any bipartite stabilizer state in 
the case of squarefree D is equivalent, up to unitaries on 
the two parts, to a collection of unentangled single-qubit 
states and maximally entangled EPR pairs. This could 
have been studied using the Schmidt decomposition, but 
the techniques used here are also needed in the following 
section. 

The central result of this paper is the tripartition 
Theorem [7] stated and proved in Section IIVI It shows 
that when D is squarefree a stabilizer state on three 
parts can be decomposed into a tensor product of single- 
qudit states, two-qudit EPR pairs and three-qudit GHZ 
states. With the help of Choi-Jamiolkowski isomorphism 
or map-state duality, this result is applied in SectionFVlto 
the stabilizer code channels where we show they can al- 
ways be decomposed into a product of a perfect quantum 
channel, a perfectly decohering channel, and a depolar- 
izing channel (not all of which need be present) . We also 
prove that the subset information groups corresponding 
to a stabilizer code channel and its complementary chan- 
nel obey a duality relation, in that one completely speci- 
fies the other. While the results are specific to stabilizer 
code channels, we show that they can also be used to pro- 
vide bounds on channel capacities for some other cases. 

Section IVTl summarizes our findings and suggests some 
directions for future research. 



II. PRELIMINARY CONCEPTS AND 
DEFINITIONS 

A. Qudit Pauli Operators 



tational basis as 
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where the addition of integers in Eq. ([T|) is modulo D. For 
a collection of n qudits we use subscripts to identify the 
corresponding Pauli operators unless otherwise stated: 
thus Zi and Xi operate on the space of qudit i. The 
Hilbert spaces of individual qudits are denoted by H^, 
and that of n qudits by H®" -.^ Hi ® 'H2 ® ■ ■ ■ ® Hn- 
Operators of the form 
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will be referred to as Pauli products, where A :— e^^^/'^'^^) 
(so A^ — to) and 7 is an integer in ^,20, the ring of integers 
modulo 2D. For a fixed n, the collection of all possible 
Pauli products in Eq. Q forms a group under operator 
multiplication, the Pauli group Vn- 

For every p G Vn, P^ is either / or — /. The order of 
a Pauli product p E Vn is defined as the smallest integer 
1 ^ a ^ D such that p" ex /. Our definition of order 
is nonstandard in that we only require the power of the 
Pauli products to be proportional to the identity. Note 
that the order of any Pauli product must divide D. 

While Vn is not abelian, it has the property that any 
two elements commute up to a phase: piP2 — '-^°'^^P2Pi, 
with ai2 an integer in Z/^ that depends on pi and p2- One 
can find subgroups of Vn that are abelian, for example 
the set of Pauli products with only powers of Z on every 
qudit. 

Proposition 1. Let A be set of mutually commuting 
Pauli products in Vn (for example, abelian subgroups of 
Vn)- Then A can have at most D" linearly independent 
elements. 

Proof. The elements of A can be viewed as D"-hy-D" 
matrices. Then it is impossible to simultaneously diago- 
nalize Z?" -I- 1 or more mutually commuting and linearly 
independent Z)"-by-£'" matrices. D 

The collection of Z?^" Pauh products in Eq. ^ with 
7 = 0, i.e. a pre- factor of 1, forms an orthonormal basis 
of £(7^*^"), the space of linear operators on H^", with 
respect to the Hilbert-Schmidt inner product 
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Most of the following preliminary concepts have been 
introduced in [I^l, [lal and we present them here again 
for completeness. We generalize the notion of Pauli oper- 
ators to higher dimensional Hilbert spaces where D ^ 2. 
The X and Z Pauli operators are defined in the compu- 



B. Single-Qudit and Two-Qudit Clifford Operators 

Having defined Pauli operators, we now generalize 
other single-qubit and two-qubit operators to D ^ 2. 



Pauli operator 
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TABLE I: The result of conjugation of Pauli operators by one- 
qudit gates F, S'°' and W . {a is the multiplicative inverse of 
a mod D and A = e^'^'/'^^'.) 



The qudit generalization of the Hadamard gate is the 
Fourier gate 
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For an invertible integer a € Zd (i-e. integer for which 
there exists a G I^d such that aa = 1 mod D), we define 
a multiplicative gate 
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(6) 



The requirement that a be invertible ensures that S^"' is 
unitary. (For D — 2, the only invertible integer is a = 1, 
hence S^"^ is just the identity.) 
Next we define the phase gate as 
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J:^^' X-iU+^)\j}{j\ if I? is even 
^J^j-;A-^(^+i)|j)(j| if i^ is odd. 
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where A = e^'^'/'^^-^^. The phase gate was first studied by 
Nielsen et al. in [1^] for the general D case. 

The three single-qudit operators defined above as well 
as the Pauli operators defined in Eq. (HI) are examples of 
Clifford unitaries, by which we mean unitaries that map 
Pauli products to Pauli products under conjugation. For 
instance, FZF'< = X and FXF^i = Z~^. The resuhs of 
conjugating the Pauli operators by F, S*^"^ and W are 
summarized in Table HI 

The generalizations to Z? ^ 2 of CP and CNOT gates 
are the Clifford unitaries 



CPi 



and 
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where qudit 1 is the control while qudit 2 is the target. 
The CP and CNOT gates are related by a local Fourier 
gate defined in Eq. (O, similar to the D = 2 case, 



Proposition 2. For D prime, let p € Vn be a Pauli 
product on n qudits [Eq. ^] and assume that p is not 
the identity on qudit 1, i.e. xi j^ or zi ^ or both. 
Then there exists a Clifford unitary U such that UpW oc 
X1I2 ■ ■ ■ In- Further if p^ — I , then it is possible to have 

UpU^ =Xih---In- 

Proof. If xi — 0, then conjugate p by the Fourier gate, 
F so that xi ^ 0. Then transform X^^ Z^" to Xf^ by 
conjugating it with the W gate a sufficient number of 
times. Next use the S''"-' gate to produce Xi. See Table 
Ufor the result of these conjugations. (Note that we relied 
on the fact that Z^ is a field when D is prime in the last 
two operations. The more general result for arbitrary D 
is studied in [lal-) If P is now the identity on all the other 
qudits i = 2, 3, . . . , n we are done. Otherwise, for each 
non-identity qudit, set the Pauli operator to X employing 
the procedure above. If at this point p is not identity on 
qudit 2, i.e. p = X1X2 • • • , then X2 can be changed to I2 
by performing CNOT12. This is repeated where needed 
so that the Pauli product is the identity on every qudit 
except qudit 1 , which proves UpW is proportional to Xi . 
If now p^ = /, any remaining phase is necessarily some 
power of w, and can be removed by conjugation with 
powers of Zi . D 



C. Stabilizer Codes and States. Partitions 

Let S C Vn be an abelian subgroup consisting of lin- 
early independent Pauli products. Then \S\ must divide 
D" and s^ ^ I for all s E S. (For prime D, every ele- 
ment except the identity is necessarily of order D so \S\ is 
always a power of D.) Given S, define the set of states, 
C := {|^) e -H®" : s|V') = \ip), ^ s e S}. It is easy 
to check that C forms a linear space, which we call the 
stabilizer code, with S its stabilizer group [32| . 

In [l^ it was shown that C and S are dual in the sense 
that one completely specifies the other and they satisfy 
the relation \S\ x dim(C) — Z?"; dim(C) is the dimension 
of C In quantum error correction literature, if dim(C) — 
D^ for some integer ^ A: ^ n, then it is customary to 
write C = [[n, fc]]£i because one can think of encoding k 
qudits in the ZJ'^-dimensional subspace contained in the 
space of n carrier qudits. Let [/ be a Clifford unitary 
and S a stabilizer group with C being its corresponding 
stabilizer code. Then S' := USW = {UsW : s e 5} is 
also a stabilizer group stabilizing the code C — {U\^) : 
IV') S C}. For a detailed review on Clifford unitaries and 
stabilizer states for arbitrary D see ^^ . 

If \S\ ~ D", then S stabilizes a unique state and we 
call it the stabilizer state, denoted by \S). The projector 
onto the state can be written as a sum of elements in S, 
as shown in \\M 
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Two simple examples of stabilizer states are the EPR 
pair and the GHZ state, expressed below for any D ^ 2 
with their respective stabilizer groups, 
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and 



D-l 



|GHZ)i23 = ^T. I*)i 



V2|i;3, 



S — (X1X2X3, Z1Z2 , Z\Z-^ 



(13) 



where the angular brackets denote the group generated 
by products of the elements in the list. 

When D is prime, S can always be generated by n 
suitably chosen group elements, S = (si, S2, . . . , Sn) such 
that the order of each Si is D. For non-prime D, one 
might need more than n generators is some cases. We call 
these group elements stabilizer generators or generators. 
Note that the set of generators is not unique - there 
are many distinct choices of generators that generate the 
same group, for example S — (siS2, S2, • • ■ , Sn)- 

Proposition 3. Let S be a stabilizer group with D" el- 
ements where D is prime. Let T = (ii,i2j ■ • ■ j^m) &e a 
subgroup ofS with D"^ elements where 1 ^ m < n. Then 
there exists a set of n ~ m elements, {tm+i, ■ ■ ■ ,tn} d S 

such that S = (^1,^2,- • • ,im,im+l,- ■ • ,tn). 

Proof. First pick an element of S not in T and call it 
tm+i- Since D is prime, the order of tm+i must be D. 
Then the set {i'C+il *' G {ti,t2, ■ ■ ■ ,t^),a G Zd},= 
{ti,t2, . . ■ ,t„i,t,n+i) is a subgroup of 5 with _D™+i el- 
ements. Repeat this incremental addition of generators 
until the set of generators generates S. D 

A stabilizer state \S) € "H®" naturally "lives" in a 
tensor product space of n qudits but one can imagine a 
coarser-grained partitioning where the n qudits are di- 
vided into two parts, labeled A and B, which we will 
call a bipartition. One can regard any state on the total 
Hilbert space as an entangled state on Ha (X" H-b and we 
will also refer to such a state as a bipartition. Triparti- 
tions and generalizations to higher number of partitions 
can be analogously defined. (Obviously partitions can 
be defined on any multi-partite state, not just stabilizer 
states.) 

A useful expression for reduced density operators of 
multipartite stabilizer states is the following. For a bi- 
partite stabilizer state \S) £ Ha ^ Hb let 
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be the elements of S equal to the identity on Hb ■ They 
form a subgroup of S (e.g., 113), and in light of Eq. ([TTj). 
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If we square both sides we see that p^ = {\Sa\/D^*)pAi 
which means that the reduced density operator of a sta- 
bilizer state has identical positive eigenvalues, so it is 
proportional to a projector. Additionally it satisfies 



rank(pA) 



D"" 



\Sa\ 



(16) 



Therefore, pA is proportional to the identity if and only 
if the subgroup Sa has only the identity element. 

Finally let \S) and \T) be stabilizer states on dis- 
tinct sets of n and m qudits with stabilizer group S — 
(si, . . . , Sn) and T = (ii, . . . , i™). Then clearly the ten- 
sor product of these states |V) = \S) ® \T) is also a 
stabilizer state with the stabilizer group 



V = (si,...,S„) (g) {ti,...,tjrt) 

= (si (g) /,..., s„ (g) /,/ (g) ii 



,/(gi™). (17) 



Conversely, given a stabilizer state |V) e "H®" (gH®™, if 
the stabilizer group can be written as a tensor product 
of two stabilizer groups, V = (si, . . . , s„) g) (ti, . . . , tm), 
then |V) = |5) g) |T), since \S) and \T) are uniquely 
determined by their respective stabilizer groups. 



D. Decomposition of Stabilizer States of 
Composite Dimensions 

Let the integer D have the prime decomposition 



d^p1^pI^---pI 



(18) 



where the pi are distinct primes and the e^ positive in- 
tegers. The following theorem is useful when the qudit 
dimension D is composite. 

Theorem 4 (Chinese Remainder Decomposition of Sta- 
bilizer State). Let A be an abelian group of linearly in- 
dependent Pauli products on n qudits, each of dimension 
D, and let Eq. (1181) be the prime decomposition of D. 

Then A is unitarily equivalent to a tensor product of 
m abelian groups in the sense that 



{U®---® U)A{U ( 
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where U is a unitary acting on the D-dimensional space 
of a single qudit, each Ai is an abelian group of linearly 
independent Pauli products on n qudits of dimension p')^ , 
and 



1^1 = |A|- 1-421- --lA 



(20) 
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The proof is in Appendix A. As stabilizer groups are 
examples of such abelian groups, one has: 

Corollary 5. Let \S) be a stabilizer state on n qudits of 
dimension D, prime decomposition given by Eq. (J18p . 



Then there exists a single- qudit unitary U such that 
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U' 



)U\S) = 



(21) 



where each \Si) is a stabilizer state on n qudits of dimen- 
sion pI^ . 

Proof. Let S denote the stabilizer group of \S) which has 
£)" hnearly independent Pauli products. Then by the 
theorem above S, an abehan group of hnearly indepen- 
dent Pauli products, is equivalent up to local unitaries 
to a tensor product of m stabilizer groups of dimensions 
pI^ , P2^ , • ■ • , Pm" I each stabilizing its own stabilizer state 

\S^). D 

When applied to an arbitrary stabilizer state \Se) on 
qudits of I? == 6, this corollary states that it is equivalent 
up to local single-qudit unitaries to a tensor product of 
two stabilizer states |iS2) {S^), one on n qubits and the 
other on n qutrits. Essentially each D = 6 qudit has 
an internal tensor product structure that can be decom- 
posed to a qubit and a qutrit. Therefore in studies of 
entanglement of stabilizer states of D = 6, it is sufficient 
to just consider qubit and qutrit stabilizer states. 

In this paper, the corollary above is used to extend 
various results on stabilizer states that hold for prime D 
to the case where D is squarefree, meaning that e^ = 1 in 
p8)) for every i. 



E. Graph States 



Let Tij = Tji be the adjacency matrix of an undirected 
graph G on n vertices with no loops {Tu =0)). Each Tjj, 
the weight of the edge connecting vertices i and j, can 
take any value in Zjj, with (as usual) Tij = in the 
absence of an edge. 

The graph state \G) is a state on n qudits of dimension 
D defined as 



|G) - n n cp 



1=1 j=i+l 



where 
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and the two-qudit gate CP ij is defined in Eq. ([S]). Note 
that the CP gates all commute with each other so there 
is no need to specify the order in which they act on the 
kets. For non-prime Z?, there are alternative ways to 
define graph states; see, e.g., |19| . 

All graph states are stabilizer states but the converse is 
not true; however it was shown in [20[ and |19j for prime 
D (and therefore also for squarefree D by Corollary [S]) 
that all stabilizer states are equivalent up to local Clifford 



unitaries to graph states. It is often more convenient to 
work with the stabilizer group (denoted by Sq) rather 
than the ket itself. For a given graph state \G) with 
adjacency matrix F, there is a canonical set of n stabilizer 
generators, {gi} given by 



■■=x,. 




for i = l,2,...,n (24) 



which of course satisfies gi\G) — \G) for all gi, so we 
have Sq = (51 , 52 , • ■ • , Qn) ■ These operators are called 
correlation operators in [21| . 



III. BIPARTITION OF QUDIT STABILIZER 

STATES 



Entanglement across bipartitions of stabilizer states 
has been studied in Section 3 of |22| and fis!]. Here we 
shall extend their result to all squarefree D ^ 2 with the 
theorem below. The entanglement of a bipartite state can 
always be studied in terms of its Schmidt decomposition 
but we present an alternative approach here because it is 
helpful in explaining the techniques that will be used in 
the Tripartition Theorem in Section IIVI 

Before stating the Bipartition Theorem let us study 
some simple stabilizer states to understand how unen- 
tangled subsystems in each part can obscure the actual 
amount of entanglement present. We shall consider the 
two stabilizer states on three qudits below and ask how 
much entanglement is present across the A-B bipartition: 



S^ — {Zai^Bi iXAiXbi,Xa2 



^ — i^Ai ^Bi ' ^Ai 'Z'A-i ^Bi 1 ■2'ai X^ 



(25) 



where the subscripts Ai,A2 denote qudits in part A and 
analogously Bi in part B. 

First observe that iS^^-* can be factorized as 
{Zj^^Zg^,XAiXBi) <8) {XA2) (see discussion at the end 
of Subsection IIICp and the unentangled qudit A2 is ir- 
relevant as far as entanglement between A and B is 
concerned. From here it is straightforward to see that 
{Za^Zbi, Xa^Xbi) stabilizes the EPR pair described in 

Eq. dni). 

In the second case, it is harder to tell how entangled the 
state is by just looking at the stabilizer group 5'-^-' , even 
though it differs only by a local unitary on part A from 
the previous state, I^S^^)) = CPaiAsI'?'^^). This tells us 
there must be some hidden unentangled subsystem in 
part A that upon removal will result in a simpler two- 
qudit stabilizer state, just like the first example. 

A systematic way to "detect" the presence of unentan- 
gled subsystems in stabilizer states is by inspecting the 
reduced density operator on each part or equivalently the 
subgroups 5a, 5b fsee Subsection HTCt . The Bipartition 
Theorem below is essentially just a formal statement that 
once all the unentangled subsystems are removed, all that 
remains is a collection of EPR pairs. 



Theorem 6 (Bipartition of stabilizer state). For square- 
free D, let \S) he a stabilizer state on n ^ 2 qudits. For 
any bipartition of \S) G Ha ® 'Hb, there exists Clifford 
unitaries Ua,Ub on each part such that UaUb\S) is a 
collection of maximally entangled EPR pairs and unen- 
tangled single- qudit states, i.e. 



UaUb\S)ab 



I EPR) ^'"Ai 



I AB 



<E> |+)t"^ <E> \+)T 



(26) 



Note that niA, niB or tuab can be zero. 



Proof. By invoking Corollary [S] we can decompose \S) 
into several stabilizer states where each of them is on n 
qudits of prime dimension. Therefore it is sufficient to 
prove the theorem only for prime D. 

If the subgroups Sa and Sb both contain only the iden- 
tity element, then by Eq. (|15|) both pA and pB are pro- 
portional to the identity. This is equivalent to S not 
containing any element that is non-trivial only in one 
part, such as Xa^Ib- This also means \S) is maximally 
entangled and therefore is equivalent to a collection of 
EPR pairs. 

Otherwise assume that Sa has at least one element, s G 
Vn not equal to the identiy, a Pauli product which acts 
non-trivially on at least one qudit in part A. Without loss 
of generality we can assume that qudit is Ai. Then by 
Proposition [2] we know there exists a Clifford operation, 
Ua such that Uji^sU\ = XaiIa2 ' ' ' Ia„ ■ 

Next consider the new stabilizer group for Ua\S) and 
choose Xai as one of the generators so that Uj^SU\ = 
{Xai , S2t ■ ■ , s„), which is always possible as shown in 
Proposition [31 Since S2 must commute with Xai, there 
cannot be any Za-^ operator in it and hence must be 
of the form S2 — X^ Pa\Ai j where Pa\Ai is some 
Pauli product on qudits A2, . . . , An^. If a — 0, then 
do nothing. Otherwise replace S2 with 53 :— X^"s2 = 
I Ax ®Pa\Aii so that the new generator is the identity on 
qudit Ai. This replacement does not change the group 
being generated. 

Repeat this procedure for all the other generators 
S3, . . . , s„. In doing so we now have a new set of gen- 
erators such that there is only one generator that is 
non-trivial on qudit Ai while all other generators have 
identity on Ai. The end result is a stabilizer group 
that can be written as a tensor product, UaSUa ~ 

{Xa,)®{s'2,...,s'J. 

Following the discussion at the end of Subsection III CI 
we can write Ua\S) = \+)ai ^ \S') where |-|-) is defined 
in Eq. (^5]) and \S') is stabilized by {s'2, . . . , s'j). In other 
words, we have extracted an unentangled subsystem from 
part A and are left with a stabilizer state with n — 1 
qudits. Repeat this process on both parts until both Sa 
and Sb contain only the identity element. This concludes 
the proof. D 

In fact this extraction of unentangled subsystems 
works for stabilizer states with any number of parts since 
we can always view the part of interest as A and all the 
other parts as B. 




Ua,Ub,Uc 



A 




FIG. 1: Example of equivalence of a tripartite graph state to 
a collection of single-qudit states, EPR pairs and GHZ state. 



IV. TRIPARTITION OF QUDIT STABILIZER 

STATES 

The problem of tripartition of qubit {D — 2) stabi- 
lizer states has been studied by Bravyi et al in 12] . They 
proved that such states are always equivalent up to uni- 
taries on each part to a collection of GHZ states, maxi- 
mally entangled EPR pairs and unentangled single-qubit 
states; see Figure [T] for a simple illustration. In the same 
paper, they also provided partial solutions to the gen- 
eral problem with more than three parties. Here we ex- 
tend their tripartition result to squarefree D ^ 2 using a 
method mentioned but not used in their paper. 

Theorem 7 (Tripartition of stabilizer state). For 
squarefree D, let \S) be a stabilizer state with n > 3 qu- 
dits. For any tripartition of \S) £ Ha ® Hb <8 He, there 
exists Clifford unitaries Ua , Ub , Uc on each part such 
that UaUbUc\S) is a collection of GHZ states, maxi- 
mally entangled EPR pairs and unentangled single-qudit 
states, i.e. 

UaUbUc\S) = |GHZ)®™-^«<^ (g) lEPR)®™-*^ 
(g) |EPR)^'"«c ^ |EPR)®"-*c 



1+)' 



\(gim_E 



(27) 



where mA,rnB,rnc,rnABi'mBC,'m-AC,'m-ABC o,re non- 
negative integers that can be zero. 

Proof. It is sufficient to prove the theorem only for prime 
D because Corollary \E\ extends the proof to squarefree 
D. The proof can be divided into three major steps. 
In Step 1 we simply repeat the unentangled subsystem 
extraction procedure explained in the previous section; 
the details will not be repeated here. As a consequence 
the three reduced density operators pA, Pb and pc are all 
proportional to the identity. In Step 2, we extract EPR 
pairs using local Clifford unitaries acting on two parts at 
a time. Thus for parts A and B, we try to find Ua, Ub 
such that UaUbIc\S) = |EPR)ab ® \S'). The identical 
procedure can be used to extract EPR pairs from A and 
C, and from B and C, so we only need to discuss how 
it works for A and B. Finally in Step 3 we prove that 
the state remaining after all these extractions must be a 
collection of GHZ states. 



We begin Step 2 with the observation that the reduced 
density operator pj^s for the combined parts A and B is 
given by 



\Sab\ 

PAB = Trc{|5)(5|} = D-("-^+"-) Y. ' 

\Sab\ 



(i) 
AB 



j^-iuA+UB) Y^ gW ^ „W 



M ^ '^'S 



•Sr I 



(28) 



,(0 J») 



are 



where Sab '■= {s G S : Trc{s} 7^ 0} and 
Pauh products on A, B respectively. There is some phase 
ambiguity in the final expression because a phase on s^ 
can be moved to s^ or vice versa, but the following proof 
does not depend on how the phase is assigned. 

STEP 2 A - If the collection {s^^} defined in Eq. ^ 

contains two elements sj^ ,s\ that do not commute, 
then at least one EPR pair can be extracted from parts 
A and B as we now show. We can always assume those 
two elements satisfy the commutation relation 



^A ^A 






(29) 



because if the phase picked up is instead some higher 
power of w, we can replace s . with an appropriate power 



Jfe) 



(3) 



oi Sj^ , since D is prime. Now transform sj^^ so that it 
is the identity on every qudit except Ai and Bi, i.e. 



,0) 
^ab 



"Ai 



"Bi' 



(30) 



by applying Proposition[2]twice, first on part A and then 
on part B, using Clifford unitaries local to each part. The 
commutation relation between s^\ s^ in Eq. ([^^ is left 
unchanged because it is invariant under conjugations by 
unitaries. 

STEP 2B - The fact that [s^]^, s^^^] = together with 
Eqs. (gH) and ^ imply that 



^AB — ^A- Z 



AiPa\Ai ® ^Bi^Bi'iBXBi 



(31) 



where p, q are some Pauli products on the remaining qu- 
dits on part A, B respectively and a, /3 G Zu- To see this, 
first observe that Eqs. (P^)) and (PH)) tell us the exponent 
on Xai is necessarily 1 while imposing no conditions on 
the exponent of Zai nor the Pauli operators on the other 

qudits in part A. Next, s]^g commuting with s\^ fixes 
the exponent of Xb^ to be 1, resulting in Eq. ([51]) . 
In Eq. (l5Tt the Z operators on qudits Ai and Bi can 

(k) 

be removed by conjugating s^^ with the phase gate W 
applied to these qudits (see Table |l|. Lastly, the Pauli 
products p and q can be set to / using methods out- 
lined in the proof of Proposition [5J The conjugations by 
CNOT gates described there with Ai , Bi being the con- 
trol qudits do not modify s^g at all. At the end of all 
these transformations we have 



^AB — ^M 



'Xf 



(32) 



STEP 2C - Keeping track of all the Clifford unitaries 
that have acted so far on the two parts and combin- 
ing them as Ua , Ub , we can write the new stabilizer 
group_as {UaUbIc)S{UaUbIcV = {h = Xa.Xb.M = 
Zj^^Zq ,^3,...) by Proposition [3l Since tj, must com- 
mute with t\ and ^2, it must be of the form 



-/3, 



-^Ai^AiPa\Ai ® ^Bi^Bi 1b\Bi 



(33) 



for some a, /3 € Z^, where p, q, r are some Pauli products 
on subsystem A less Ai, B less Bi and C respectively, 
reusing arguments that produced Eq. (pij) . Next replace 
the generator t^ by fg — t^ t^^t^ so that t'^ has identities 
on qudits Ai,Bi. The same can be done for all other 
generators, ^4, . . . , i„. 

STEP 2D - Given this, we can write the new stabilizer 
group as {Xa^Xb^.Z^^Z]^^) (g) (ig, ■ • ■ ,i^): an EPR pair 
stabilized by {Xa^Xbi, Zj^ ^b^) and a stabilizer state 
\S') stabilized by (tg, . . . ,i^)- The reduced density op- 
erators p'a and p'b corresponding to \S') are again pro- 
portional to the identity operators on the parts of A and 
B that remain after the extraction. Were it otherwise 
in the case of p^, the reduced density operator on the 
(full) system A corresponding to the state stabilized by 
{XAiXBi,Zj^^Zg^) (g) (ig,...,^^) would not be propor- 
tional to the identity, contradicting the fact that after 
Step 1 (extraction of unentangled subsystems) , and thus 
at the beginning of Step 2, pA was proportional to the 
identity; of course the same applies to p'^. 

STEP 2E - Next examine the expansion of p^^ putting 

\S') in (|28p . If some of the s]^' do not commute with each 
other a further extraction is possible, and one can repeat 
the process until all of these operators commute. Now 
apply the same extraction process to A and C, and then 
to B and C, until all EPR pairs have been extracted. 

The final step is showing that the tripartite stabilizer 
state after extracting all unentangled states and EPR 
pairs is a collection of GHZ states. 

STEP 3A - We prove that after the preceding extrac- 
tions have been carried out, the three parts must have 
the same number of qudits, ua = nB = nc- Since pc is 
proportional to the identity, \S) is maximally entangled 
across the C-AB cut. Using its Schmidt form allows the 
the projector on \S) to be written as 



\S){S\ = 



1 



£)"£• 



Y. V){J\c®\<i^i){(kAAB, 



(34) 



with {|0i)} a set of orthonormal kets in Ha ® Hb- The 
{|z)(j|c} are a set of D^"^ linearly independent opera- 
tors, as are the {\4>i){(f)j\AB}, which means the operator 
Schmidt rank of \S){S\ is D^"^. 

From Eq. pTI) the projector can also be expressed as 
a sum of the Z?" linearly independent stabilizer elements 



8 



in S 



fe=i 

= —T 



„(fe) 



^C^) ^ .C') 



(35) 



fc=i 



ik) (k) 

where the r^^ and r^ are PauH products on AB and 
C, respectively. In general, the collection {rj^ } is not 
linearly independent; e.g., each r^*"'' that belongs to the 
subgroup Sab satisfies r|-, oc Ic- Indeed, two elements 
r^'^^ and r^'^ belong to the same coset of iS^b if and only if 



„(fc) 



oc r, 



(0 



Therefore the number of cosets of Sab is the 



number of linearly independent elements in the collection 

Vc }k^=i- ^^ shall now prove that this number is Z?^"'^. 
We can re-express Eq. (j35p as a sum over linearly in- 
dependent Pauli products on part C 



D-'^c 



|5)(5|(x ^ 



(fe') (*;') 

1c '^ ^AB ' 



(36) 



„(fe') 



where {xab } ^''^ sums of Pauli products on parts A and 
B. We know there must be D^""^ hnearly independent 
terms since that is the operator Schmidt rank and thus 
none of the Xab can vanish. Hence the number of cosets 
of Sab is D^"^ . Thus by Lagrange's theorem, 

\Sab\ - |5|/i?2"C ^ Jjn^+ns-nc^ (37) 

We now go on a digression to show that the collection 
of Pauli products on subsystem A, {sa}, on the right 
side of Eq. 



is linearly independent when pB oc I b', 
similarly, the collection {s^ } is linearly independent if 
Pa oc Ia- Because they are Pauli products it suffices to 
show that no two of them are proportional in order to 
demonstrate linear independence. Assume the contrary, 
that Sa cx: Sa for some j ^ k. Since the {s54;g} are 

group elements, Sab must have an inverse, I s^g ) , 
which inserted in Eq. (P5)l yields 



^AB i ^AB 



s^P®s^i>] {s'^'^s's') 



oc I A 






(38) 



The final term cannot be proportional to the identity, as 

(?) (k) 

that would imply s^g oc s\g. But these are elements 
of S, so they must be linearly independent for j ^ k. 
Therefore Sb contains an element that is not the identity, 
contradicting the fact that pB oc Ib- Hence it cannot be 
the case that Sa oc Sa for j ^ k. 

Thus the collection {s^ } contains |iSab| linearly in- 
dependent and mutually commuting elements, and by 
Proposition [T] this means that \Sab\ ^ D"^ . The same 



argument applies to the collection {sg }, so \Sab\ ^ 
Z?"^. Combining these inequalities with Eq. ([571) . it fol- 
lows that both nA and n^ cannot be larger than nc- 
Identical arguments applied to different pairs of subsys- 
tems implies that 



nA = n.B = nc, 



(39) 



and, using Eq. ^, \Sab\ = \Sbc\ = \Sac\ = D''^ . 

STEP 3B - If Sbc contains only the identity element, 
then ua = nB = nc = and no GHZ state can be ex- 
tracted. Otherwise ua ^ 1 and Sbc has at least one 
non-trivial element which we will label as ti. By Propo- 
sition [5] it can be transformed to 



U =Z^ 



Z? 



(40) 



Now consider the group Sab- We showed previ- 
ously that the operators {sg } defined in Eq. (^5)) are 
both linearly independent and mutually commuting, and 
[{•^B }l ~ -D"^. Hence, if there is an element p € Vub 
that commutes with every element in {sg }, it must be- 
long to {s^ } up to a phase. This is because, by Proposi- 
tion [1] sets of mutually commuting Pauli products on ns 
qudits cannot have more than D"^ linearly independent 
elements. 

Note that f i commuting with every element of Sab (as 
they are all just elements of S) implies Zg commutes 

with every s^ since ti is identity on part A and every 
s G Sab is identity on part C . Then by the argument in 
the previous paragraph, there exists an element in Sab 
of the form pA ® Zb^ ■ Let ^2 denote this element and 
note that pA 7^ Ia as otherwise this would contradict the 
assumption that pB is proportional to the identity. Then 
by Proposition [21 there exists unitary transformations 
such that 



h 



z 



Ai 



<z 



Bi 



(41) 



and they do not affect ti since all the operations are done 
on part A. 

STEP 3C - Since each of the D^'^c p^uii products 

for part C must appear in Eq. (|36p — none of Xab can 
vanish — there exists an element in S which satisfies 

ik') 

Qq oc XciIc\Ci on the C subsystem. We call that ele- 
ment t^. The most general form it can have, given that 
it has to commute with ti = Zq Z'^ and ^2 = Z~^ Z^ , 

is 



P 



t3 = Xa,Z'X^pa\a, «) Xb.Z'^^qbxb, <E) Xc, 



(42) 



where p, q are Pauli products on the remaining qudits 
in part A, B respectively and a,/? G "Ed', see STEP 2B 
for the explanation. Finally we transform this element 
to is = Xa-iXbiXci without modifying ii,t2 by using 
the techniques described in STEP 2B and in the proof of 
Proposition [21 

STEP 3D - Invoking Proposition [3] and let- 
ting Ua,Ub,Uc denote all the unitary oper- 
ations we have made so far, we can write 



iUAUBUc)SiUAUBUc)^ = (il = Zg^Zc^,t2 = 

Z^^Zg^,t3 ^ XAiXBj^Xci,t4,.. . ,tn)- As was done 
in STEP 2C and the proof of Bipartition Theo- 
rem, the generators t4,...,t„ can all be made to 
be the identity on qudits Ai, Bi, Ci simultane- 
ously, so we can write UaUbUcS{UaUbUc)'^ — 
{Zj,^ZcIz^IZb,,Xa,Xb,Xc,) ® {t'^,...,t'J. There- 
fore we have a GHZ state, see Eq. (fT3|) . on qudits 
yli , i?i , Ci , tensored with a state on a system which has 
one fewer qudit in each part. This extraction process 
can be repeated until ua = nB = nc = 0. D 



APPLICATION IN QUANTUM ERROR 
CORRECTION 



In this section we apply the Tripartition Theorem to 
solve a problem in the area of quantum error correc- 
tion. It allows us to understand the structure of a class 
of quantum channels derived from qudit stabilizer codes 
(see Subsection IIICp of prime D by decomposing them 
to a tensor products of perfect quantum channels, per- 
fectly decohering channels and completely depolarizing 
channels. 

The connection between tripartite stabilizer states and 
stabilizer codes is worked out in Subsection IV Al us- 
ing Choi-Jamiolkowski isomorphism or map-state dual- 
ity. This leads to a class of channels which we call sta- 
bilizer code channels defined in Subsection IV Bl whose 
decomposition into simple channels is the topic of Sub- 
section IV CI A duality between the subset information 
groups of such a channel and its complementary channel 
is demonstrated in Subsection IV Dl 



Isomorphism Between Stabilizer Codes and 
Stabilizer States 



Let Ha and Hb be two Hilbert spaces, and {|»)a} an 
orthonormal basis for T-La- Then there is a one-to-one 
correspondence between a linear map M from "Ha and 
Hb and a ket \M) on the tensor product Ha SJ Hb con- 
veniently expressed in Dirac notation as changing bras to 
kets or vice versa: 



M = ^iA)(*u, im) = ^ia: 



Wa. 



(43) 



Here the {|/3i)} are elements of Hb, in general neither or- 
thogonal nor normalized, uniquely determined by M or 
by \M) as the case may be. This Choi-Jamiolkowski iso- 
morphism 33] depends on the choice of the basis {|i)yi}; 
in what follows this will always be the computational ba- 
sis. For convenience we introduce a normalization factor 
of y/dA in defining the following isomorphism and its in- 
verse: 



with M and \M) related by Eq. (l43l) . It is straightforward 
to show that 



^UbMUa) 

<P-' [{Ub ® Ua)\M)] 



(U^UbMM), 
Ub^-'{\M))U^, 



(45) 



where the transpose T refers to the {1*}^} basis, and Ua 
and Ub are any operators on Ha and Hb , respectively. 
Next assume that dA ^ rfs and define an isometry 
V -.Ha-^Hb '-iS 



i=0 



(46) 



where {|wi)} is an orthonormal basis spanning a dA- 
dimensional subspace of Hb that we call the coding space. 
The isomorphism carries it into 



\V):^HV) 



1 



d-A 
1=0 



N)a, 



(47) 



where, since V is an isometry, pA — Tvb {\V){V\} = 
lAJdA- One can think of \V), where Ha and Hb enter 
on an equal footing, as an atemporal representation of 
the isometry |23J . 

The following lemma uses the Choi-Jamiolkowski iso- 
morphism to relate isometrics corresponding to stabilizer 
codes to stabilizer states. 

Lemma 8 (Isomorphism between stabilizer code and bi- 
partite stabilizer state). Let D be the dimension of any 
one of the qudits. The two statements below are true for 
any prime D 



1. Given an [[n,k]]D stabilizer code (Subsection \II C\ ) 
that defines an isometry V : H'^^ -> "H®", the iso- 
morphism $ in Eq. (j44[) carries V to a stabilizer 
state \V) on k -\- n qudits. 



$(M) = |M)/VdI, $"'(|M» 



'dA M, 



2. Let \S)ab be a stabilizer state on k-\-n qudits, where 
the first k qudits constitute part A and the remain- 
ing n qudits constitute part B, and assume that 
Tt::b{\S){S\] — Ia/L)^- Then the inverse isomor- 
phism ^~^ in Eq. (1441) carries \S)ab to an isometry 
V : Ha -^ Hb whose image or coding space is an 
[[n,k]]r) stabilizer code. 

Proof. The proof can be simplified by noting that, as 
shown in [l^ [20|, any stabilizer code when D is prime 
is equivalent to an additive graph code, up to products 
of single-qudit Clifford unitaries. In turn it was shown 
in [15| and [13| that an additive [[n, A;]]d graph code 
when D is prime can be described using a graph state 
\G) e ?^®" and an abehan group J" = {qJ^o"^ with 
D'' linearly independent Pauli products in Vn, composed 
only of Z operators. The coding space is spanned by 
the set {ci|G)} of mutually orthogonal kets. Recall from 
Subsection III CI that any graph state \G) on n qudits can 
(44) be fully specified by its stabilizer group, (gi, . . . , gn) with 
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D" elements, see Eq. ((24)) . Similarly, the coding group 
T can be generated by fc suitably chosen group elements, 
^ = (/i: /2, ■ • ■ , /fe), so the coding space is spanned by 
{/i' ■■■fk" |G)} for zi, i2, . . . , ifc = 0, 1, . . . , i? - 1. 
To prove statement 1, define the isometry 



D-l 



£1-1 



^ = E---E/"---/^iG)^(*i" 



■«fcU 



(48) 



ii=0 



ifc=0 



that maps "Ka = 'H'^'^, a collection of fc qudits, into the 
coding space, and 



\y) = 



1 



D-l 



^ 41=0 



_D-1 

■E 

U.=0 



l«l • • • jfc)A «) /l' 



(49) 



the corresponding ket on 'Ka ® "Hb as per Eq. p4l) . 

We shall show that \V) is a stabilizer state by exhibit- 
ing the k + n stabilizer generators. The first n are derived 
from the generators {gj}, for j — 1, 2, . . . , n of the stabi- 
lizer group of |G) , now regarded as operators on 'Ha^'Hb, 
so that 



ilA(^9j)\V) 



1 



D-l 



rnk /-^ 

^ ii=0 



D-1 

■E 

ifc=0 



il 



_. D-1 £)-l 
^ -^ ii=0 ife=0 



■ik). 



■Ik) a 



■/fc1G)B, (50) 



where the phases result from commuting the gj with the 
/; operators: 

9jfi = ^'^" fi9] ioi j ^l,...,n and I ^l,...,k; (51) 

with the f3ji integers in Z^. The phases can be removed 
using an appropriate Pauli product of Z operators: 



'z,^''---z: 



-Pik 



9j)\V)^\V) forj-l,...,n. 



(52) 



That these n generators mutually commute follows from 
the fact that {gj} is a mutually commuting set. 
The remaining k stabilizer generators are simply 

^A, ®(/r')B for/ = l,...,fc. (53) 

They obviously commute among themselves. It is not 
hard to show that they leave \V) in Eq. (|49)) invariant, 
and commute with the n previous generators on the left 
side of Eq. ([5^ . Since we have constructed k + n linearly 
independent commuting Pauli operators that leave \V) 
invariant, and thus generate a stabilizer group, \V) is a 
stabilizer state. 

To prove statement 2 of the lemma we use the fact that 
for prime Z?, any stabilizer state \S) is equivalent up to 



single-qudit Clifford unitaries, thus a choice of basis for 
the individual qudits, to a graph state. Hence without 
loss of generality we can assume that |iS) is a graph state 
on "Ha ® T~iB written in the form, see Eq. 



\s)ab = n ^^' 



i+>r- 



(54) 



i<j 



Let CP^ be the product of those CP gates on the right 
side that act only on qudits in A, and let \G)b — 
CPs|-l-)f" be the graph state on "Hs resulting from the 
CP gates that act only on qudits in _B, denoted by CP^. 
The action of the remaining CP gates that connect A 
and B qudits can be written out explicitly in terms of Z 
operators, see Eq. (|5]), to obtain 



\S) 



k n+A; 

cp^ I n n CP 

.i=i]=k+i 



1 



D-l 










'\G)i 



■Ik) A 



■fk'\G)i 



where 



n-\-k 



fi 



n <'^ 



for ? = 1, 2, . . . , fc. 



(55) 



(56) 



]=k+i 
and so 

n+k D-1 

n cp[/-en')(*'I®//' 
j=k+i 



for I = 1, 



il=0 



(57) 

Comparing the last line of Eq. ([55]) with Eq. (|i^ yields 
D-l u-i 

$-M(CP^)t|5)] = ^...^/r.../^1G)B 

n=0 ik=0 

(ii---JfcU, (58) 

which is an isometry corresponding to an additive graph 
code with graph state \G)b & 7^®" and coding group 
C = (/i, . . . , /fc). Therefore, using Eq. ([iS]) 

'i>-'(|5)) = ^---^C/s/r---/^1G)B 

ii=0 ife=0 

{ti---ik\A[{CPA)Y, (59) 

where the transpose is taken in the computational ba- 
sis. The final factor is simply a unitary transforma- 
tion on the input and thus does not change its image, 
which is the graph or stabilizer code spanned by the 
/r • ■ ■ fk\G)B- That these last are a collection of D*-' 
mutually orthogonal kets follows from the assumption 
that Tri3{|5)(5|} = Ia/D^, and the fact that the fi- 
nal equality in Eq. (|55p is a Schmidt decomposition of 
\S). U 
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B. Stabilizer Code Channels and Subset 
Information Groups 



of simple channels, which has important implications for 
the properties of Qb and Qc- 



Consider an isometry V : I-La — "H 



qj®k 



n 



(8in 



corre- 



sponding to an [[n, /c]]_d stabilizer code where D is prime 
and the n output qudits are partitioned into two dis- 
joint non-empty subsets, B and C. Such bipartitions of 
stabilizer/graph codes have been studied in JJJ and ^4] . 
(Qubit stabilizer codes with input qudits partitioned into 
two parts have also been studied in [25] but in this paper 
we will only consider bipartitions of the output qudits.) 
We can think of the B qudits as the output of a direct 
quantum channel, and the C qudits either as the environ- 
ment or as the output of complementary channel, with 
corresponding superoperators 

£b{p) :- Tvc{VpV^} Ec{p) := TrB{VpV^}. (60) 

We shall refer to channels derived in this way from stabi- 
lizer codes as stabilizer code channels. The analysis be- 
low applies to any stabilizer code regardless of its error 
correction properties such as code distance. 

In Section V of ^^, we studied the information- 
carrying capacity of stabilizer code channels and to this 
end introduced the subset information group 



Gb ■■={peVk \£b{p)^0}, 



(61) 



a subgroup of the Pauli group on the k input qudits. It 
was shown that Qb is a group, and its elements satisfy 
the isomorphism 



£b{p) £B{q) = c£B{pq) yp,q^GB, 



(62) 



where c is an appropriately chosen positive constant in- 
dependent of p and q. 

We also presented an efficient algorithm to find Gb 
given the isometry V (defined by an additive graph code) 
and the subset B, by solving a set of linear equations 
modulo D. The subset information group Gc for the 
complementary channel can be defined in the same way. 
The previous work discussed additive graph codes, but 
the results apply to stabilizer codes as well since they are 
equivalent up to local unitaries. 

One can think of the group Gb or the operator alge- 
bra that it spans, a subalgebra of the algebra C{Ha) of 
operators on the channel input, as representing the in- 
formation that is perfectly transmitted from the input to 
the output T-Lb of the channel by £b ■ Since it is present in 
the output, this information can be perfectly recovered, 
which is to say mapped to a Hilbert space "H^ isomor- 
phic to Ha , by a recovery operation (recovery channel) 
TZ : C{'Hb) -^ CCH'a), as shown in [ij]. The recovery 
operation has the property that [TZ o Eb]'' (x) — x, where 
X is any Pauli product in Gb or any operator in the sub- 
algebra that it spans. Thus the last is an example of a 
correctable algebra, as defined in [14| in their formalism 
of operator algebra quantum error correction. 

In the following subsections we show that stabilizer 
code channels can be decomposed into tensor products 



C. Tensor Product Structure of Stabilizer Code 
Channels 

Before stating the main result in Theorem [9l let us 
indicate by means of some simple examples its main idea, 
which is that the isomorphism proven in Lemma [5] with 
Theorem [7] imply stabilizer code channels have a very 
simple structure. 

First consider the stabilizer state, \V) = |EPR)ai_Bi, 
Eq. ([H]), which by Eqs. gS]) and ^ corresponds to the 
isometry 



D-l 



'^-\\V}) = Vepr^Y. 



\VB^{^\A^ 



(63) 



j=0 



for a perfect quantum channel from Hai to Hbi , with 
quantum (and classical) channel capacity equal to log2 D. 
The subset information group contains every Pauli oper- 
ator on qudit Ai, i.e. Gb — {^Iai,Xai, Za^}- 

Next consider the tripartite state \V) = \EPR)aiBi ® 
|EPR)b2C'i • Tracing out part C, Ea. ([5(H) . yields the chan- 
nel £b with the same Gb = {^Iai,Xa-i, Za-^) as be- 
fore, while the complementary channel £c is the com- 
pletely noisy or completely depolarizing channel whose 
subset information group is (multiples of) the identity, 
Gc = {^Iai)- Therefore, attaching the state |EPR)b2Ci 
to that of the previous example increases the dimen- 
sion of the output Hilbert space while leaving the Ha 
to Hb channel unchanged. This is not surprising, since 
the |EPR)b2C'j^ part has nothing to do with the input 
Hilbert space Ha- 

As a third example, the tripartite |GHZ) state from 

Eq. inD 

1 ^"' 
\V) = \GRZ) A,B,c, = ^^ E l*)^i l*)^! I*)-^! (64) 



is carried by the inverse map in Eq. (j44p to the isometry 
'^~\\V)) = l^GHZ = E N)bJ*)c,(*U,. (65) 

i=0 

In this case, £b and £c are perfectly decohering channels 
whose Kraus representation is 

£b{p)^£c{p)^J2\')('\p\')('\^ (66) 

i=0 

generalizing the qubit phase-flip channel (see Chap. 8 
of [55) to arbitrary D. The quantum channel capacity is 
zero while the classical channel capacity is log2 D. The 
subset information groups for both channels are identical: 
Gb = Gc = (A/ai,^Ai)- 
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B 



C 



\J 



FIG. 2: Decomposition of the isometry of a [[7, 3]]i3 stabilizer 
code into simple isometrics as proven in Theorem |9l The 
dotted line in the diagram on the right represents an EPR 
pair between parts B and C. 



The following theorem states that the isometry of any 
stabilizer code with a bipartition defined on the output 
qudits is equivalent to a tensor product of isometrics of 
the form Vepr and Vqhz ■ See Figure [5] for an example 
of a decomposition of an isometry of a [[7, 3]]_d stabilizer 
code with ub — 3 and nc = 4. 

Theorem 9 (Tensor product structure of stabilizer code 
isometrics). For prime D, let V : Ha -^ Hb (Xi He be 
an isometry corresponding to an [[n,k]]D stabilizer code 
with a B-C bipartition of the n output qudits. Then up to 
unitaries Ua, Ub, Uc on Ha, Hb and He, V is a tensor 
product, 



UbUcVUa 



rr^niABC 

"GHZ 

® |EPR)« 



K 



EPR 



rr^rriAc 
*^EPR 



+ 



(67) 



where rriABC, "mAB, "itt- AC, 'mBC "^Bi itt-C ^'^e non- 
negative integers, and Vepr and Vqhz are defined in 
Eqs. (pl and ({6^1). 

Proof. Use Lemma |8] to map the stabilizer code isometry 
y to a stabilizer state ^{V) = \V) on k + n qudits. Use 
Theorem [7] to express the result, up to local unitaries, 
in the form given in Eq. (P7)) . Apply to this the inverse 
map ^~^, noting that unitaries can be pulled outside, as 
shown in Eq. (gS]) . D 

An immediate corollary is that stabilizer code chan- 
nels, up to unitaries on the input and output spaces, 
are always tensor products of just three types of sim- 
ple channel; (i) Perfect quantum channel, (ii) Perfectly 
decohering or phase-flip channel, (iii) Completely noisy 
or completely depolarizing channel. This determines the 
quantum and classical capacities since, e.g., the quantum 
capacity of £b is just the number of EPR pairs that can 
be extracted from the A-B bipartition of \V) multiplied 
by log2 D, while the classical capacity is the quantum ca- 
pacity plus the number of GHZ states linking A, B, and 
C. 

Hence it follows that if £ is a stabilizer code chan- 
nel and its classical and quantum channel capacities are 
known, then one can deduce the decomposition of the 
isometry as in Eq. (j67p and Figure O since the dimension 
of the input is known. From here it is straightforward to 



work out the classical and quantum channel capacities of 
the complementary channel. 



D. Duality of Subset Information Groups, Qb and 

Gc 

Here we shall demonstrate how the tensor product 
structure of stabilizer code channels enables us to eas- 
ily determine the subset information groups, Gb and Qc. 
We also prove a duality relation between the Gb and Qc 
which implies that each of them is determined by the 
other. 

Recall that the subset information groups Gb and Gc 
are subgroups of the Pauli group Vk on the k input qu- 
dits, and that the centralizer Centg(y^) of any subset A 
of elements of a group G consists of all elements of G that 
commute with every member of A. 

Theorem 10. For D prime, let V : Ha ~^ Hb ® He be 
an isometry corresponding to an [[n, /c]]/) stabilizer code, 
with Vk the Pauli group on the input qudits. Assume a 
B-C bipartition of the n output qudits is defined, with Gb 
and Gc the corresponding subset information groups as 
defined in Eq. (j6ip and its counterpart with B replaced 
with C , for the stabilizer code channels £b and £c ■ 
Then Gb and Gc ol^R dual in the sense that 



Gb ^Geni-p^^{Gc) and Gc ^ Cent-p^iGB) 



(68) 



i.e., each is the centralizer of the other in the group Vk- 

Proof. First we show that the duality relation holds for 
Fepr and for Vghz, as defined in Eqs. (p5|) and (I55|) re- 
spectively. For V^pR*^' , we have Gb = (A/^i , Xa^ , Za^) 
from the previous subsection. The complementary chan- 
nel is Scip) = Tr b{VpV^ = Tr{p}, and by Eq. (EH) we 
have Gc — {^Iai)- Thus Gb is the full Pauh group, 
and Gc consists of multiples of the identity, so they 
are centralizers of each other in the Pauli group on 
qudit Ai. Next, for the isometry Vqj^^ we have 

Gb = Gc ^ {^IaiiZai) and again the duality is satis- 
fied. 

But Theorem [S] says that V is equivalent up to local 
unitaries to tensor products of these two types of elemen- 
tary isometrics, where we ignore the single-qudits states 
on i3, C and the EPR pairs between parts B and C as 
they are not involved in the transmission of infornration. 
It is straightforward to show that since the Pauli group 
for part A and the information groups are themselves 
tensor products, the latter are again centralizers of each 
other. Since the duality is invariant under conjugation 
by unitaries, the local unitaries play no role. D 



Here is an illustrative example. Suppose that 



V = V. 



EPR 



v; 



EPR 



J/As—^B2C2 

^GHZ 



(69) 
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Then the subset information groups have the same tensor 
products structure as the isometrics above 

Gb - (A/) (E) {Xa, ,Za,)® {I a,) <E} (Za,) 

= (A/, Xai ,Zai, Za3 ) , 
Gc = {XI) ® {I A,) ® {Xa,,Za,) ^ (Za,) 

= {XI,Xa,,Za,,Za,). (70) 

While the results presented in this section only hold 
for stabilizer code isometrics, we now show that they 
can sometimes be used to derive bounds for more gen- 
eral isometrics V, with channels £b and £c defined as 
in Eq. (j60p for some bipartition of the output qudits. 
Given a general coding space V, one can ask if there 
are stabilizer subspaces or subcodes contained in V. If 
these subcodes are not one-dimensional subspaces, then 
meaningful lower bounds on channel capacities can be 
calculated for £b and £c- 

For example consider the nonadditive (or non- 
stabilizer) qubit graph code denoted by ((5, 6, 2))2, men- 
tioned in [l5| and [2g . (It was first described in \2r\ , but 
not using the graph code formalism.) Also assume the 
following bipartition of the five output qubits, B = {1,2} 
and C = {3,4,5}. The six-dimensional coding space is 
spanned by 

V = spB.n{\G),ZiZ2Zi\G),Z2Z3Z5\G), 

ZiZ^ZA\G),Z2ZiZ5\G),Z^Z^Z^\G)}, (71) 

where \G) is a five-qubit graph code stabilized by 
the group (X1Z2Z5, Z1X2Z2,, Z2X3Z4, ZsX^Z^ZiZ^X^) 
and the corresponding graph is a pentagon 15]. 

Obviously, V contains the two-dimensional stabilizer 
subcode Vq C V, 



Vo = span{\G),ZiZ2Zi\G)}. 



(72) 



By applying the techniques used in proving Lemma |8] and 
Theorem [H] one can show that Vo corresponds to a one 
qubit quantum channel from A to G, which means that 
the quantum channel capacity of £c is greater or equal 
to log2(2). Similarly, by considering the subcode 



Vi=span{|G),Z2Z3Z5|G)}, 



(73) 



one can identify a GHZ state in the corresponding tripar- 
tition, and thus deduce that the classical channel capacity 
of £b and of £c is at least logj (2) . 

Working from the other direction, one can also consider 
stabilizer codes that contain V as a subcode. Provided 
these stabilizer codes are not the whole Hilbert space, 
useful upper bounds on channel capacities of £3 and £c 
can be calculated. 



these allows stabilizer states for composite D to be ex- 
pressed as tensor products of stabilizer states associated 
with the different prime factors of D. This is a valuable 
technical tool used later in the paper, but also a helpful 
conceptual tool as it allows more complicated cases to 
be "pulled apart" into simpler situations. For example, 
a graph state for D = 6 can be regarded as the tensor 
product of D = 2 and D = 3 graph states. 

Our main result. Theorem [3 that a tripartite stabi- 
lizer state can be considered the tensor product of single- 
qudit states, two-qudit EPR pairs and three-qudit GHZ 
states, generalizes the D = 2 (qubit) result in j^ to any 
squarefree D > 2. This allows us to provide a very sim- 
ple and essentially complete characterization, TheoremlHl 
up to unitaries on the input and output, of channels con- 
structed from stabilizer quantum codes. Knowledge of 
the information carrying properties of such a channel 
leads immediately tells one the properties of the com- 
plementary channel, and there is a simple relationship 
between the corresponding subset information groups in- 
troduced in [13|. In some cases one can use these results 
to put bounds on capacities of other types of channel. 

There are various directions in which one might hope 
to extend these results. We have encountered techni- 
cal difficulties in attempting to generalize our triparti- 
tion theorem from squarefree to arbitrary composite D, 
where it may no longer be true. It follows from Corollary 
[5] that it is sufficient to resolve the situation in which D 
is a prime power, so if that could be solved one could 
have results that apply for general D. Can any of our 
results be extended beyond the narrow confines of sta- 
bilizer states? The fact that the proofs depend heavily 
on group-theoretical properties makes this seem unlikely, 
but it would certainly be of interest to understand what 
it is that makes stabilizer states so special, and stabilizer 
quantum codes so useful. 

Another possible direction is to move on from tripar- 
titions to those involving four or more parts. Here the 
results in [3] for qubits suggest a situation that is dis- 
tinctly more complicated than found for bipartitions and 
tripartitions, and Z? > 2 is unlikely to be simpler. But it 
would still merit study. 
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VI. CONCLUSION 

The most important results of our paper are those in 
Corollary [SI Theorem [71 and Theorem [5) The first of 



Appendix A: Proof of Theorem [4] 

Proof. It is crucial to first note that there are two dif- 
ferent notions of tensor product used in the statement 
of the theorem and in this proof. The first and more 
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obvious notion is the tensor product of n qudits, each 
of dimension D. The second notion is defined on the 
space of each qudit of dimension D, which is isomorphic 
to a tensor product of m spaces on qudits of dimensions 
di,d2, . ■ ■ , dm ■ 

First of ah, we show that Pauli operators on a single 
qudit of dimension D are equivalent to tensor products 
of m Pauli products on the constituent qudits. We will 
be relying primarily on the Chinese Remainder Theorem 
which states that there exists a ring isomorphism between 
Zjj and the tensor product of rings, Z^j (8) Z^a (Xi • • • €5 Z^^ 
whenever D has the prime decomposition in Eq. (J18p . 

For brevity it is sufficient to prove the theorem for the 



case oi D = did2 where di — pi^ and d2 



P2 



■Pn 



be- 



cause d2 can subsequently be further decomposed by in- 
duction. The Chinese Remainder ring isomorphism map, 

(/) : Zd — > Zdj (g) Zd2 is defined as 



0(a) = {(t>i{a),(t>2ia)) 



(Al) 



where 



0i(a) :— a mod di. (A2) 

The inverse map ip^^ : Z^^ (g) Z^a -^ Z/j is given by 

(/i^ {ai,a2) '■= airid2 + a2r2di mod D (A3) 

with Ti := [D /di)~^ mod di being constants that depend 
only on di,d2 and not inputs ai,a2- Note that r^ is al- 
ways coprime to di for i = 1,2. 

Next we show the mapping (f> induces a unitary trans- 
formation from a Hilbert space of dimension £> to a tensor 
product space of qudits of dimensions di , ^2 . We define 
the action of the unitary U : Hd -^ 'Hdi ® T~l-d2 on the 
D-dimensional basis kets as 

l(\a) = |0i(a)) 102(a)) for a = 0, 1, . . . , L> - 1. 

(A4) 



with 0i's defined in Eq. (|A2I) . The fact that (p is bijective 
guarantees that {|0i(a)) (g) |(/>2(a))}£^o^ spans the space 
■Hdi '^'Hd2- 

The result of conjugating Pauli operator X by this 
unitary is 

UXU^ = Y, \Ma)){Ma + 1)1 ® \Ma)}{Ma + 1)| 

0=0 
D-1 

= E \Ma)){Ma) + 1| \Ma)){Ma) + 1| 

a=0 
di-ld2-l 

= 51 H |ai)(ai + 1| <8) \a2){a2 + 1| = ATi ® X2. 

ai—0 ao—0 

(A5) 

In the last line we simply replaced the sum over all ele- 
ments of Z^) with the sum over all elements in Z^^ (8 Zj;^ . 
Note the subscript on X here denotes two difi^erent qu- 
dits, each with different dimension. 



Next the Z operator is transformed as 

D-l 



UZU^ - Y, wl0i(a))(0i(a)| (i> |02(a))(02(a)| 

a=0 
di-1 

= J^u;*"('^i''^^)|ai)(ai|®|a2)(a2| 
) 

di~l 

Y e'"'"'"'/*|a.)(a,| = Z[i ® Zl,\ (A6) 



ai=0 
2 di-1 



i=l ai=0 



where r^'s are defined in Eq. (jA3|) . 

We are now ready to prove the theorem. Let A be 
generated by fc elements, A = {g^^\ . . . , g^'^^) and the 
generators can always be chosen such that 



nordcr(5W) = |^|, 



(A7) 



where the order of every gi must be a divisor of D, see 
Subsection III Al for our nonstandard definition of order. 
The requirement that y^ be a collection of linearly inde- 
pendent Pauli products implies every generator satisfies 

9^ '"'=/. 

Consider an arbitrary generator, g^"^' and define 5 = 
order((7^'^). Let 5 — P1P2 • • 'P^ be the prime decom- 
position of 5 with the same piS as Eq. (|18l) . and set 

5i — Pi , 62 = P2'"Pm\ so (5 = &i52- Then 82 is 
coprime to 5i and also to pi. Next define the unitary 
\] = U®- ■ -^U from Eq. (IA4[) acting on all the n qudits. 
By Eqs. (|A5p and (|A6p . conjugating the first generator 
with U produces 



Ug(^)U^ = qi 



'92 



(A8) 



where gi , 52 are Pauli products on n qudits of dimension 
di , d2 respectively. Next we claim that 



order((7i) = 5i and order((72) 



(A9) 



To prove this, recall that (g('))''i'^2 — / which implies 
qi ^ oc /i. Since 82 is coprime to 5i and also to di, the 
unique multiplicative inverse of $2 mod di exists. Then 

it follows that qi "^ '^ — qi (x Ii. Therefore, the order 
of qi must be a divisor of 5i and by the same reasoning, 
the order of 92 rnust be a divisor of 82. The orders cannot 
be less than 5i , 82 respectively as that would imply the 
order of g^'-' is less than 6. 

Having proven Eq. (|A9p . we define Pauli products of 
qudits of dimension di and ^2, hi and /ij respectively 
as 









(AlO) 
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where Hi :— {S/5i) ^ mod S.^. Next for all ai € T^Si and 
a2 G Z^j , we can rewrite Eq. (jAlOp as 



(/.«)" «(4'T^-K^ut) 



Qi^ii52+a2A»2<5i 



(All) 



Observe that the exponent on the right side is just the 
inverse Chinese Remainder map of ai and a2 inEq. (|A3p 
applied to this situation, so 



Ug«Ut) = (/if>)®(4*>) 



(A12) 



Decomposing every generator into two generators on 
subsystems of dimensions di and d2 gives us 



Uy4U^ = (Ug(i)u^...,U<?('='U^) 



.(1) 



,W\ 



,(1) 



.W\ 



A (8) ^2- (A13) 



That Ai 



{h\ 



(1) 



.,/i^''^) and A2 



{h^\ 



,W\ 



form collections of mutually commuting Pauli products 
are consequences of the fact that {5^*^} are mutually 
commuting and that conjugation by U does not change 
the commutation relation. Finally, Eq. (|A13[) also tells 
us there are as many linearly independent elements in 
A as there are in Ai ® A2, which implies that |^| — 
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